SPHERICALLY SYMMETRIC SOLUTIONS OF A BOUNDARY VALUE 

PROBLEM FOR MONOPOLES 



ANTONELLA MARINI AND LORENZO SADUN 

Abstract. In this paper we study spherically symmetric monopoles, which are critical points 
for the Yang-Mills-Higgs functional over a disk in 3 dimensions, with prescribed degree and 
covariant constant at the boundary. This is a 3-dimensional gauge-theory generalization of the 
Ginzburg-Landau model in 2 dimensions. 

1. Introduction 

In this paper we treat a 3-dimensional analogue of the vortex equations in 2 dimensions 
and look for solutions with spherical symmetry as described in §2. The domain considered is 
a 3-dimensional disk and we prescribe the degree of the monopole at the boundary. Unlike 
in 2 dimensions, where Ginzburg-Landau-type functionals appear either with or without gauge 
potentials, the problem in 3 dimensions is well-posed only in gauge theory. Without a curvature 
term in the action, a minimizing sequence for the action would yield a trivial limit [Q]. In the 
presence of a gauge potential, this problem is well-posed, natural and has physical meaning. 

The most general Yang-Mills-Higgs functional takes the form 

(1.1) yMn{A,4) = '-\\F\\l, + ^mnh + ^11 |<^P - a'wl^ , 

for appropriate constants e, p, A and a. Working on M^, one usually applies a rescaling of (p, a 
rescaling of space, and a rescaling of the action to set e = p = a = 1, so the action functional 
depends on a single parameter, A. On the unit ball, however, we cannot rescale space, so we 
can only eliminate two of the four parameters. We set p = a = 1, and obtain a 2-paramater 
family of functionals 

(1.2) yMn,,x{A, 4>) = |||F||i.(53) + ^WDAHh^BS) + ^11 \cP\' - l||i2(B3) . 

(Alternatively, we could work on a sphere of radius R. One can then rescale to set e = 1, at 
the cost of varying R. We then obtain a 2-parameter family of functionals indexed by A and 
R.) We know from the general theory for monopoles (cf. ^ for e = 1, A > 0) that there exists 
a minimum for this functional which satisfies the Euler Lagrange equations 

e*DA*F= [Da(P, 4>] 

(1-3) A , 

*DA*DA(p=-m^ -1)4 
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and suitable boundary conditions on dB'^ = (cf. Section 4 and and is smooth. In 

this paper we prove the existence, and describe the form, of sphericahy symmetric solutions to 
these equations. 

We note that, even for A = 0, these are not solutions to the Bogomolnyi equations found in 
[^. The Bogomolnyi solutions are obtained only in the limit A — > 0, i? ^ oo, or equivalently 
A ^ 0, e ^ 0. 



2. Spherically symmetric connections, monopoles, and gauge transformations 

We work on the trivial principal S'C/(2)-bundle P = x SU{2) and its associated vector 
bundles. ^ is a connection on P, which can be viewed as a 1-form on B"^ with values in 
su{2). The Higgs field (p is a section of the adjoint bundle, i.e., a map (/) : B'^ ^ su{2). Here 
su(2) = {X e M2X2 ■■ trX = 0; X + = 0} isthe Lie algebra of SUi2). We identify su{2) 
with the imaginary quaternions ImM. = {xii + X2j + x^k : {xi,X2, X3) G M'^} by identifying the 
matrices j) , ;) , with fc, respectively, and extending this mapping to 

a Lie-algebra isomorphism. It is also convenient to define the 3-vector of Lie-algebra elements 

= {i,j,k). 

The symmetry group S0{3) acts on pairs {A, (p), simultaneously rotating the 3-dimensional 
base space and the 3-dimensional Lie algebra. That is, the triple k) transform in the same 
way as the triple {xi, X2, X3), so quantities such as a ■ x are invariant. 

We are interested in finding Yang-Mills-Higgs fields {A, 0) which are invariant under this 
group action. To this purpose, one needs to specify the value of the connection one-form 
A : B"^ A^(B'^) ImM. and of the Higgs field (j) : B^ ImM. at one point of each group 
orbit (on the base) and impose invariance under the isotropy group of that point. We find it 
convenient to fix the values of {A, (p) on the slice L = {x G B"^ : X2 = X3 = 0, xi G (0, 1]}. The 
isotropy group at (3;,0,0) G L is S0{2), i.e. rotations about the xi-axis (and about the i axis 
in the Lie-algebra). For the Higgs field (p one has in general 

(j){x, 0, 0) = (pi{x)i + (p2{x)j + ipz{x)k G ImM ~ su{2) , 

where (pi{x) , I = 1,2,3 are real-valued functions. Imposing invariance under S0{2) forces 
^2ix) = ^-iix) = for all x G (0, 1]. Applying the action of SO (3), one obtains the symmetric 
form of the Higgs field (p 

(2.1) (I) = ■ x = ^ {xii + X2j + xsk) , 

with r = \x\, X = {xi,X2,X3) G M'^, a = k) G ImM. 



Yang-Mills-Higgs, symmetric bundles 3 

An su(2)-valued connection on the slice L is given in general by 

A{x, 0, 0) = aii{x)i dxi + ai2{x)i dx2 + ai3(x)i dx^ + a2i{x)j dxi + 

(2.2) 

+ a22{x)kdx2 + a23{x)j dxs + a3i{x)kdxi + a32{x)kdx2 + +a33{x)k dxs 

Imposing S0(2)-invariance yields 012 = ais = 021 = 031 = 0, 022 = 033, and 023 = —032- Thus 
the final form of the 50(2)— invariant connection evaluated at the points (x,0, 0) G L is 

A{x, 0, 0) = a{x)i dxi + b{x){j dx2 + k dx^) + c{x){k dx2 — j dx^) . 

Transporting this slice via the SO (3) group action on one obtains the invariant version 

(2.3) A{xi,X2, X3) = a ■ dx H ^ dx) • (x x <?) H —a ■ {x x dx) , 

where"x" denotes the cross product of vectors and dx = {dxi,dx2,dx3). 

At this point there is still some gauge freedom available to further specify the connection 
A, namely that provided by "symmetric" gauge transformations. Such a transformation g is 
determined by its values on the slice L, which must be S'0(2)-invariant. This yields i g{x, 0, 0) = 
g{x, 0, 0) i, thus g{x, 0, 0) = exp {i h[x)). Therefore, symmetric gauge transformations are of the 
type 

a ■ X, 



(2.4) g{xi,X2,X'i) = exp{h{r) 

where h{r) is an arbitrary function of the radius r. Performing such a gauge transformation 
does not change the form (l2.lD of the Higgs field (j). However, A transforms nontrivially. In 
particular, setting h{r) = j ^^j^ dr, exactly cancels the a-piece in g~^dg + g^^Ag. Thus, one 
can impose a{r) = 0. The final version of A is then 

(2.5) A(xi,X2,X3) = — :t-(x X dx) ■ (x X a) -\ 5-cT • (x X dx) . 

The only gauge freedom remaining is from the constant of integration in the indefinite integral 
J f^M (cf. Section 4). 

Note that the f3 and 7 terms have opposite parities. The isometry x — > —x sends /3 to — /? 
but sends 7 to +7. 

3. The Yang-Mills-Higgs functional on spherigally symmetric configurations 

In this section we explicity compute the Yang-Mills-Higgs functional, and the resulting equa- 
tions of motion, for symmetric pairs {A, (p). Our connection is a sum of two terms, 

A = B + C, 
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where 



(3.1) ^ 

and 



P{r) j3{r) 
B = — ^ ^2;) • (x X cj) = — 3~[('' — — xix^k — xiX2j]dxi + 



P(r) f3(r) 

+ ~::3~[('^^ ~ ^2)i - X2X3k - XiX2i]dX2 H —[{1^'^ - xl)k - XiXsi - X2X3j]dX3 , 



C = "^^J a ■ (x X dx) = '^^ J (x^j — X2k)dxi + 

(3.2) ^ ^ 

t(^) t(^) 

+ —^{xik - X3i)dx2 + —^{x2i - xij)dx3 
Preliminary computations 
The curvature of A is given by 

F = dA + AAA = dB + dC + BAB + CAC, 

since the cross term BaC + CABis identically zero. 
Computing F12, the various nonzero terms are 

P' 

(3.3) {dB)i2 = -^{-X2i + xij) 

(3.4) (dC-)n = + + ^i)) 

(3.5) [B A B)i2 = [xix^i + X2X3j + xp] 

(3.6) (C A C)i2 = [xixsi + X2X3j + xjk] . 

By rotational symmetry, the contributions to F13 and F23 are similar. 
The covariant derivative of the monopole 4) = ^^^a • x is given by 

(3.7) DA<t) = d(t>+[B,ct)] + [C,(t)], 
where 

(3.8) d(t)='^d^ + ^d{^) , 



(3.9) [B, (p] = B(f) — (pB = -^(— xsj + X2k) dxi + cyclic permutations, 

2'yip 2 2 

(3.10) [C, = C(j) — (pC = — 3~[('^ — — X1X3A; — xiX2i]rfxi + cyclic permutations. 

At this point we are ready to compute the three terms involved in the Yang-Mills-Higgs 
functional (1.2). They are 

1) = F A*F 

2) \D(p\'^ = D(j)A*D(j) 
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3) m'-lf. 

Computation of 1): One shows easily that 

dB A *dC = dC A *dB = , 
dB A *{B AB) = {B AB)A *dB = , 
dB A *{C AC) = {C AC) A*dB = 0, 

thus 

\Ff =\dB\'^ + + dC A *{B AB) + dC A *{C A C) + {B A B) A *{B A B) 

+ {BAB)A *{C A C) + (C7 A C) A *{C A C) 

^2 ^4 ^4 ^4 ^4 

?'2 o^/2 



2/3'^ 27'^ 4(/?^ + 7^ + 7)^ 



^2 ij-*^ 

Computation of 2): One shows that 

fi(?i A 0] = (/>] A = , 
[B, (/.] A 0] = [C, 0] A 0] = , 

thus 

(3.12) = |(i</.|2 + |[5,(/,]|2 + |[C,(/)]|2 + 2#A*[C,(^] 

(3.13) = f/V24V8^ + 8^ + 8^- 
Computation of 3): One easily obtains 

Collecting terms, the Yang-Mills-Higgs functional calculated on spherically symmetric con- 
figurations is given by 

(3 14) SyMHij, ^) j\ 2e + + ^(/3' + + 7)') + ^ V + 

4. Further gauge transformations and the Euler-Lagrange equations 

We want to search for absolute minima of the functional ( 3.14| ) among all finite-action spher- 
ically symmetric configurations (/3, 7,(^7) . The following theorem restricts the possibilities: 

Theorem 4.1. If the functional ( 3.14| ) has a minimum, then this minimum is achieved by 
functions {p,^,<p) with [3 identically zero and 7(0) = 0. 
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Proof. First we find a gauge transformations that yields /3(0) = 7(0) = 0. In that gauge, we 
then show that minimization requires to be constant, hence for /3 to be identically zero. 

If SyMTi.{(3,^,ip) < 00, then (3 and 7 much approach well-defined limits as r — > 0, and 
/?^(0) + 7^(0) + 7(0) = 0. If /?(0) and 7(0) are not already zero, we let 6 be the argument of 
the complex number /3(0) + ^7(0), and define 



g = exp(6 

\x\ ' 



(4.1) 
Then 

(4.2) dg(xx,X2,x-i) = sin(6') 

and evaluating on our slice gives 
(4.3) 



xxi^X23^x-ik ■ (a\^'^ 

expw — j = cos(y j + sm(C7) . 

\x\ \x\ 



a ■ dx , , 1 > ^ ^ 
— r^i h di-r^jO' • X 

\x\ \x\ 



dg{r,0,0) = — —{jdx2 + kdx^). 



Our transformed connection on the slice is then 
(4.4) 

[g-^dg + g-^Ag\{r,Q,Q) 



cos(26') (3{r) ^ sin(26')7(r) ^ cos{9) sin(6l) 



sin(26l)^(r) cos(26l)7(r) sin2(6l) 



+ 



/3(r) 'y('") 

I (j dx2 + k dxs) -\ (kdx2 — j dx^)- 

r r 



{j dx2 + kdx'i) 
{kdx2 — j dx^) 



Plugging in the values of sin(6'), cos{9), etc., gives 

(/32(0) - 7'(0))/?(r) + 2/?(0)7(0)7(r) + /3(0)7(0) 



(4.5) 
(4.6) 



/3(r) 
7(r) 



/32(0)+72(0) 
(-2/3(0)7(0))/?(r) + (/32(0) - 7'(0))7(r) - 7^(0) 



/32(0) + 72(0) 

As r — > 0, both these terms go to zero, since /3^(0) + 7^(0) + 7(0) = 0. 

Having set /3(0) = 7(0) = 0, we now show that /3 must be identically zero. We choose polar 
coordinates in the (/3, 7 + ^)— plane: 

{(3 = 1/ cos t 
1 . , 
7 + - = 1/ sm t 



In these coordinates, the functional (3.14) becomes 



(4.8) 



SyMn{i^,t) =4tt [ [2e 
Jo 

8(^V^ + Ar2(v92 



/2 , ,/2 2 



1)2 ]dr , 
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with i^^(O) = J. The only dependence on t is in the t'^z^^ term, which is minimized by setting 
t = constant. But then cot(t) = -tt niust be constant, and equal to f.^^ i = 0, so /3 is 
identically zero. □ 

We may therefore restrict our attention to the functional 

(4.9) 5(7, ^) ='^^ j\ 2e + ^(t' + t)') + ^V" + V(l + 27)' + Ar2((^2 _ ^^^2 ^ 

The Euler-Lagrange equations for this functional together with the appropriate boundary con- 
ditions are then 

on (0, 1) 

\2 o\ .2 



(4.10) 



7 



V(l + 27)-4(7' + 7)(l + 27) 



7(1) = 4 



on (0, 1) 



7(0) = 

W{l)=+l (or(^(l) = -l) 
The boundary conditions above come directly from the variational principle. In fact, to cancel 
the boundary terms one needs to either restrict the space of connections to those with prescribed 
boundary data, or to impose Neumann-type boundary conditions. The boundary condition 
(p{l) = ±1 comes from |(/)| = 1 on dB^ and 7(1) = — ^ comes from 

' a ■ X 



(4.11) 



(Z)0), = (#), + [A„ , 



^1 + 2j)(l) d 







on dB^, where the subscript r denotes tangential components (cf. |3|). 

Observation: After some computation, the system ( 4.10| ) could also be obtained by impos- 
ing spherical symmetry in (1.3), thus critical symmetric points for the action ( |4.9| ) are symmetric 
critical points (not necessarily minima) for ( |1.2| ). This is also known apriori from the "principle 
of symmetric criticality" 

5. Existence of spherically symmetric monopoles 
Our basic existence result is 

Theorem 5.1. For all values o/A > 0, e > 0, there exists a symmetric solution of 

'e*DA*F=[DA(t>,<P] on B^ 



(5.1) 



A 



* Da * DacI> 

{Dct>)r = 

}ip\ = l ondB^ 



on dB^ 



1) 



on B-^ 
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Observation: These equations do not reduce to the Bogomolnyi equations, even if A = 0. 
The Bogolmolnyi argument involves an integration by parts; on a finite domain, this results in 
a boundary contribution 

Proof. Because of the derivative terms in the action, the natural space for 7 (denoted TC-^) is 
H^{0, 1), while the natural space Tl^ for 99 is the weighted Sobolev space H^{{0,l),r'^dr). By 
the Sobolev embedding theorem, functions in TC^ are continuous on [0, 1]. Functions in are 
continuous on (0,1], but may not have a limit at r = 0. We may therefore apply boundary 
conditions to 7 at r = and at r = 1, and to (/? at r = 1. 
Let 

(5.2) = {(7, ip)Gn^xn^ : 7(1) = -i, 7(0) = 0, ip{l) = 1.} 



(5.3) <^(r) 



The action functional (4.9) is well-defined on J^, and is finite whenever ip is bounded. In 
particular, fi = InfjrS is finite. We follow the direct method in the calculus of variations. 
That is, take a minimizing sequence for 5, show that it converges weakly in and then show 
that the weak limit minimizes the action and so solves the Euler-Lagrange equations. 

Let (7n,'/3„) be a minimizing sequence for S. Since A > 0, the action is not increased if we 
make the replacement 

'-1, ifv?(r)<-l; 
V9(r), if - 1 < (^(r) < 1; 
1, if Lp{r) > 1. 

As a result, we can assume that each 99„(r) is bounded in magnitude by 1. Under these 
circumstances, the sequence (7n,y'n) is bounded in C Ti-y x 'H^p. However, balls in Ti^ are 
weakly compact, as are balls in Ti^p, so the pair (7^, ipn) converges weakly in TL-^ x Ti^p to a limit 

(Too, Voo)- 

By Sobolev, 7ra(r) and (/?„(?") converge pointwise to 7oo('") and ipoo{r), so the limiting values 
7(0), 7(1), and 99(1) are preserved, and (700, ^00) G ^- Moreover, terms in Si^jn, ^n) that don't 
involve derivatives converge to the corresponding terms in 5(7oo,y'oo)- The derivative terms 
are quadratic, hence weakly semicontinuous. As a result, 5(7oo,'/'oo) is bounded above by fj,, 
and therefore must equal /u. 

Showing that 700 and (foo satisfy the Euler-Lagrange equations ([4.10|) is then a standard 
exercise in the calculus of variations. Smoothness of (700,9500) away from r = follows by 



elliptic regularity of the equations ( 4.10 ). Smoothness at r = follows from regular singular- 
point analysis, combined with the fact that both functions are bounded (see §6 for details). This 
in turn implies that the connection and Higgs field {A, (p) constructed from (700, V'oo) comprise a 
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smooth, symmetric classical solution to the PDE system ( ^.ll ). (Alternatively, one can establish 
regularity of {A, (p) from the ellipticity of the PDE system ( |5.lD , since we are working in a gauge 
with d*A = 0.) 

□ 

6. Regular Singular Point Analysis 

In §5 we demonstrated the existence of symmetric monopoles for arbitrary A > and e > 0. 
In this section we explore their form near the regular singular point of the equations ( [4.1[1| ), 
namely r = 0. 



Theorem 6.1. Let {'y,(p) be a bounded finite-action solution to the ODE system (4.10) for 
some fixed e > and A > 0. Then there exist constants ai and 62 such that 

(p{r) =air + O(r^); 

(6.4) 7(r) =62r2 + 0(r^); 

7'(r) =2b2r + 0{r^y, 

near r = 0. In particular, ip{0) = 7'(0) = 0. 

Proof. We begin with the equation for 99, which we write as 

(6.5) ^" + ^-^ = ^(7 + 7^) + 2A,,(^^-l). 

Since ip is bounded and 7(0) = 0, the terms on the right hand side are less singular than 
those on the left hand side, and to leading order ip resembles the solution to the homogeneous 
equation 

(6.6) ^ +___=0. 

The general solution to this equation istp = air+a^2^~'^ ■ However, ip is bounded by assumption, 



so a_2 = 0. Thus the solution to (|6.5|) is, to leading order, air. 



Next we turn to the equation for 7, namely 

(6.7) 7" - ^ = + 27) + ^(27' + 37). 

Again, since 99 is bounded and 7(0) = 0, this may be viewed as a perturbation of the homoge- 
neous linear equation 

(6.8) 7" - 5 = 0, 

whose solution is 7 = 6_ir^^ + h2r'^ . Since 7 is bounded, 6_i = 0. Thus our solution to (|6 
is, to leading order, h^r'^ . 
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With these basic results, we can estimate the right hand sides of ( |6.5| ) and ( |6.7D . The right 
hand side of ( |6.5| ) is 0{r), which gives an 0{r) correction to 99", hence an O(r^) correction to 
(p. The right hand side of ( |6.7| ) is O(r^), thus giving an 0(r'^) correction to 7' and an 0{r^) 
correction to 7. □ 

One can do an expansion for ip and 7 in powers of r. Indeed, only odd powers contribute to 
ip and only even powers contribute to 7. This is seen by induction. By Theorem 6.1, p is odd 
and 7 is even through order r^. However, if 99 is odd and 7 is even through order r^, then the 



right hand sides of (6^) and ( |67^ ) are odd and even, respectively, to order r^, and so p and 7 



are odd and even, respectively, to order r^~^'^. Thus ip and 7 are odd and even to all orders in 



r. 



We can therefore write an asymptotic expansion 
9) " 



a„,r 



n' ! 



Plugging this expansion into equations ( |6.5D and ( |6.7D and equating coefficients of r" ^ gives 
recursion relations of the form 

{2k){2k + 3)a2fc+i = algebraic expression involving ai, 62, • • • > ^2fc, 

(6.10) 

(2A; + l){2k — 2)b2k = algebraic expression involving ai, 62, • • • > ^2fc-2- 

These relations do not constrain ai or 62- However, once we have ai and 621 the remaining 
coefficients are determined. The first few are 

03 =(4ai62 - Aai)/5; 

64 = (36i + e-^al) /5; 

05 ={4aib4 + 40362 + 40162 + A(o? - a3))/14; 
fee = (6i + 36264 + e"H«i«3 + «?^2)) /7; 

07 =[4(oi66 + 0364 + 0562 + 0362 + 2016264) + A(3aia3 - a5)]/27; 
^^'^^^ bs = [3(26^64 + bl + 26266) + e-\al + 20105 + 20^64 + 401O362)] /27; 

09 =[4(0163 + 0366 + 0564 + 0762 + 056I + 2036264 + 0164 + 2016265) 

+ A(3af 05 + 3oio| - 07)] /44; 
610 = [3(6^66 + 626I + 6263 + 6465) 

+ e~"'^(oi07 + 03O5 + O166 + 0^62 + 201O364 + 201O562)] /22. 
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7. Symmetries and Stability 

The action functional and the resulting Euler-Lagrange equations are invariant under two 
natural symmetries. 

(7.12) ip{r) -(p{r) 7(r) +-f{r); 

(7.13) (p{r) ^ ip{r) 7(r) — > — 1 — 7(r). 

The first symmetry ( |7.12| ) comes from the isometry x —x of B"^, which of course respects 
rotational symmetry. Since a ■ x is odd and a ■ {x x dx) is even, pulling the pair [A, 0) back by 
this isometry flips the sign of ip while preserving 7. Using this symmetry, we can fix the sign 
of 93(1), which we henceforth take to be positive. 

The second symmetry ( [7.13| ) is a gauge transformation by {a ■ x)/r. This is of the form of 
(4.2), with = 7r/2. From (4.4) it is clear that this transformation sends 7 to —1 — 7 without 
generating a (3 term or changing ip. Applying this to a connection with 7(0) = yields a new 
connection with 7(0) = —1. This connection has finite action but is singular at the origin, 
reflecting the singular gauge transformation that generated it. 

We now consider stability properties of the ODE system ( [4.10| ). These ODEs have several 
fixed points, namely 

(7.14) (7,(^) = (-i,l), (-i,-l), (-i,0), (0,0) or (-1,0). 

The point (— l) is related to (— ^,l) by the symmetry ([7.12|) , while (—1,0) is related to 
(0,0) by ( 7.13| ), so we do not need to study these. What remains is (— ^, l), (~^)0)i and (0,0). 
For the 7 = — 1/2 fixed points, we define (5 = 7 + 1/2, and the equation for 7 becomes 

(7.15) ' +g 



e 4r2 

The fixed point 7 = —1/2 is stable for (p = 1 when < e/4, but is unstable if > e/4. This 
defines a natural length scale to the problem, namely \/e/2. We should expect our solutions to 
behave qualitatively differently for r less than or greater than this length scale. Of course, if 
e > 4, then all radii r are less than this length scale. In the case of = 0, the value 7 = 1/2 is 
always stable, regardless of e or r. 

Near 7 = —1/2, the equation for 93 becomes 



(7.16) ^" + 2'^ = 2^(% + \{^^-l] 



The behavior of this fixed point depends on the value of (p. Near 99 = we have 
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which is stable for all positive values of r. Near ip = 1, however, we write (p = 1 + ( and have 



(rCr = 4A(rC) + 0(C') + 0(,52). 

This is unstable as long as A > 0, and has natural length scale I/VaX. 

To summarize, the fixed point (—1/2, 0) is stable, while the fixed point (—1/2, 1) is unstable. 
If < e/4, then there is one unstable mode, corresponding to growth of <^ — 1. If > e/4, 
then there are two unstable modes, one for (p and one for 7. 

Finally, there is the fixed point (0,0). Near (0,0) we have 

7" =27/r^ + higher order, 

(7.17) 

ip" + 2(p' /r =2(p{r ^ — A) + higher order. 

This fixed point is always unstable, with 7 growing rapidly, (p will grow exponentially if r < 
1/\/A, and will oscillate if r > 1/y/X. 




-0.6 1 1 1 ' 1 ' 1 ' 1 ' 1 -0.6 1 1 1 1 ' 1 ' 1 ' 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

e = 0.1 e=10 
Figure 1. Trajectories with fixed e. 



Yang-Mills-Higgs, symmetric bundles 13 

8. Numerical Results and Qualitative Analysis 

For any fixed e and A, and given ai and b2, one can in principle integrate the differential 
equations out to r = 1. In practice, numerical errors due to the discretization of the interval 
[0, 1] can be very bad near the origin, due to the singular nature of the ODE system there. A 
better method is to use the power series ( |6.9| ) in a neighborhood of the origin and to numerically 
integrate from there. In a discretization of 10, 000 points, we use the power series out to r = 0.01, 
or 100 lattice spacings from the origin. 

In this way we get a pair (7(1), (/'(I)) each (ai, 62)- Using Newton's method, we then find 
values of (01,62) such that (7(1), (/?(1)) = (—1/2,1). Table 1 lists the correct values of oi and 
62 for several values of e and A. 

The resulting functions ip{r) and 7(r) are sketched in Figures 1 and 2. Figure 1 shows the 
functions for different values of A and e fixed at 0.1 or at 10. Figure 2 is similar, only with A 
fixed and e variable. In each case the positive function is 93 and the negative function is 7. 
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Epsilon = 10 
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\ 










— ■ — ^ — 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



A = A = 30. 

Figure 2: Trajectories with A fixed. 



From these figures several qualitative features are clear. Although ip depends significantly 
on both e and A, 7 is practically independent of A, especially when e is large. The length scale 
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e 


A 


ai 


62 


0.1 





2.82909077 


-4.47460232 


0.1 


1 


3.14773551 


-4.92072556 


0.1 


3 


3.62692766 


-5.57110938 


0.1 


10 


4.62892 107 


-6.81947999 


0.1 


30 


6.19274693 


-8.46474894 


0.3 





2.01904955 


-1.88549902 


0.3 


1 


2.26118176 


-2.04994984 


0.3 


3 


2.66517994 


-2.31673622 


0.3 


10 


3.59550462 


-2.86817001 


0.3 


30 


5.12510342 


-3.58374045 


1 





1.67098122 


-1.02894746 


1 


1 


1.85973704 


-1.07504639 


1 


3 


2.19572981 


-1.15577783 


1 


10 


3.04898441 


-1.34041824 


1 


30 


4.55384341 


-1.58910470 


3 





1.57081044 


-0.80615986 


3 


1 


1.74184236 


-0.82078859 


3 


3 


2.05156143 


-0.84695210 


3 


10 


2.86750186 


-0.90924487 


3 


30 


4.35776101 


-0.99551235 


10 





1.53622287 


-0.73146686 


10 


1 


1.70099654 


-0.73576432 


10 


3 


2.00102288 


-0.74350147 


10 


10 


2.80198139 


-0.76219107 


10 


30 


4.28571713 


-0.78838676 



Table 1. Taylor coefHcients (01,62) for various values of (e, A). 



on which 7 changes from to —1/2 is the smaller of s/e and 1. The length scale on which ip 
changes from to 1 is the smallest of ^/e, 1/\/A, and 1. Thus changing A has the greatest effect 
when A is greater than 1, while changing e has the greatest effect when e < 1. 
The source code for these numerical results can be obtained from the authors. 

Acknowledgements 
The authors thank Prof. Karen K. Uhlenbeck for useful conversations on this subject. 



Yang-Mills-Higgs, symmetric bundles 



15 



References 

[1] A. Jaffe, C. Taubes, Vortices and Monopoles, Progress in Physics; 2, Birkhauser, Boston, 1980. 

[2] A. Marini, A boundary value problem for monopoles over a 3-dimensional disk, Preprint. 
[3] A. Marini, Regularity theory for the generalized Neumann problem for Yang-Mills Connections- Non- 
trivial examples in dimensions 3 and 4, Mathemathische Annalen, 317, 173-193 (2000);. 6. 
[4] R. S. Palais, The principle of symmetric criticality, Comm. Math. Phys. 69 (1979) 19-30. 

DiPARTiMENTO Di Matematica, Universita' di L'Aquila, 67100 L'Aquila, Italy 
E-mail address: marini@univaq.it 

Department of Mathematics, University of Texas, Austin, TX 
E-mail address: sadunOmath.. utexas.edu 



